Interplay between Zeeman Coupling and Swap Action in Spin-based Quantum 
Computer Models: Error Correction in Inhomogeneous Magnetic Fields 
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We consider theoretically the interplay between Zeeman coupling and exchange-induced swap action 
in spin-based quantum dot quantum computer models in the presence of inhomogeneous magnetic 
fields, which are invariably present in real systems. We estimate quantitatively swap errors caused 
by the inhomogeneous field, establishing that error correction would, in principle, be possible in the 
presence of non-uniform magnetic fields in realistic structures. 
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In recent years quantum computing has attracted 
widespread attention. The computers based on the prin- 
ciples of quantum mechanics, such as quantum paral- 
lelism and entanglement, promise to deliver results much 
faster than classical computers in certain tasks such as 
factoring and searching 0. The impetus for construct- 
ing real quantum computer architectures arose from the 
seminal results establishing that quantum error cor- 
rection is theoretically possible and therefore decoherence 
is not an insurmountable barrier as was assumed earlier. 
There have been many proposals for quantum computer 
(QC) architectures based on various physical two-level 
systems, such as those using laser-cooled trapped ions 
HQ, photons or atoms trapped in cavities nuclear 
spins in bulk solutions |7j] and crystal lattices j£j , spins of 
electrons trapped in quantum dots or donors [ pd| , 

spins of donor nuclei in Si |12|, superconducting devices 
[ p"3| |, and others. The minimal requirements for a QC 
architecture are the existence of fundamental quantum 
bits (qubit) and the ability to carry out single- and two- 
qubit operations such as swap and controlled-NOT, as 
well as suitable hardwares for reading (input) and writing 
(output or measurement) operations. Among the many 
obstacles facing the successful demonstration of nontriv- 
ial quantum computation in specific QC hardwares, the 
most daunting are the problem of quantum decoherence 
and the difficulty in achieving precise control over the 
various unitary operations necessary for quantum com- 
putation. In this Letter we consider the theoretical issue 
of controlling the swap operation in the proposed solid 
state QC architecture involving quantum dot spin en- 
tanglement, with externally-controlled electrostatic gates 
and magnetic fields providing the unitary operation. 

Initial work on QC hardwares concentrated on 
atomic/molecular systems partly because well-defined 
single-qubit operations are comparatively easier to con- 
trol in atomic systems such as trapped ions §,[§. Success- 
ful QC architectures will eventually require many qubits 
working in parallel, which would be difficult to achieve in 
atomic systems. Proposed solid state QC architectures 



are, in principle, scalable to many qubits, although no 
one has yet successfully demonstrated a controlled single- 
qubit operation in solid state QC systems. This is a curi- 
ous dichotomy in the current QC hardware research — the 
architectures with demonstrated single-qubit operations 
are difficult to scale up while the presumably scalable 
solid state architectures have not yet been able to demon- 
strate single-qubit control. It is thus crucial to explore 
the challenges facing coherent control of qubits in solid 
state structures, particularly the issue of possible error 
corrections in realistic systems. Among various micro- 
scopic degrees of freedom that have been considered for 
the role of qubits in solid state QC architectures, spins of 
electrons or nuclei are natural candidates because of their 
well-defined Hilbert spaces and their relatively long deco- 
herence time compared to the orbital degrees of freedom. 
In the proposed spin-based QC architectures, exchange 
interaction plays a fundamental role of establishing two- 
qubit entanglement J^,||,|ll],[l2| , while Zeeman coupling to 
an external magnetic field provides various single-qubit 
operations. Generally Zeeman coupling is treated sepa- 
rately from the exchange interaction because the 
latter originates from the Coulomb interaction (due to 
Pauli principle) while the former is purely a spin effect, 
and the two interaction terms in the Hamiltonian com- 
mute with each other if the Zeeman term is homogeneous. 
However, solid state heterostructures are intrinsically in- 
homogeneous and magnetic imperfections or impurities 
are likely to be present, leading to inhomogeneous stray 
magnetic fields. Furthermore, parallel pulse schemes, in 
which exchange interaction and inhomogeneous Zeeman 
coupling arc present simultaneously, have been proposed 
to expedite operations of spin-based QCs Therefore, 
it is necessary to explore the interplay between exchange 
interaction and Zeeman coupling in relation to spin-based 
QC architectures, as we do here. 

In this Letter we focus on the effects of Zeeman cou- 
pling when the external magnetic field is spatially in- 
homogeneous. Zeeman coupling is generally neglected in 
these studies Hjlg] based partly on the assumption that 
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the effects are small and do not lead to any qualitative 
changes. While it is true in uniform fields, this assump- 
tion may not hold in an inhomogeneous magnetic field. 
The key motivations for our study are the observations 
that spatially inhomogeneous fields are intrinsic to spin- 
based solid state QC architectures and that for parallel 
schemes the magnetic fields at the locations of the two 
spins are generally different |l6|, which means the total 
field must be inhomogeneous. If only the single-spin evo- 
lution is involved, effects of local magnetic field can be 
corrected in such single-qubit operations by techniques 
such as spin echoes |T^]. However, if two-spin entangled 
evolution is also involved, it is not clear whether errors 
caused by inhomogeneous fields can still be eliminated, 
i.e. whether such errors lie within the current QC error 
correction constraints. 

As an example we study how Zeeman coupling may 
affect the proposed operations of spin-based quantum 
dot quantum computers (QDQC). Our results, how- 
ever, are quite general and can be applied to other 
spin-based models with finite magnetic fields. In par- 
ticular, we study the effect of a finite inhomogeneous 
magnetic field on swap actions. For most of the spin- 
based schemes, swap action U sw , in which two spins ex- 
change their states, is one of the most basic operations. 
It is used to construct conditional phase shifts (CPS), 



U CPS = e l ~ a ^e 
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'Usw (with a being the 



Pauli matrices for spins), which can then be converted to 
controlled-NOT (CNOT) easily §. In addition, swaps 
are used to move spin states around ||[l2) so that an 
arbitrary pair of spins can be brought into controlled en- 
tanglement, which is essential to quantum computation 
jl| . It is therefore important to investigate how Zeeman 
splitting will affect the swap gate and what are the con- 
sequences if these effects are non-trivial. 



In a single envelope function approach |15|, the Hamil- 
tonian for two electrons trapped in a lateral double- 
quantum-dot is 



+g*H B B(ri) ■ Si] 
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Notice that here the magnetic field has a spatial depen- 
dence. The Hamiltonian can be expanded in a basis 
of two-spin eigenstates (spin singlet and triplet states). 
If the magnetic field is uniform, the Zeeman coupling 
depends only on the total spin along the field direc- 
tion and commutes with the total spin operator S 2 , so 
that singlet and triplet states remain eigenstates of the 
two-electron system. On the other hand, if the mag- 
netic field is inhomogeneous, it destroys the symmetry 
between the two spins, and the singlet state and one 
of the triplet states (the one with S a = 0, where n 



is the field direction) couple. For example the cou- 
pling between the lowest singlet and triplet states (in the 
Heitler-London approximation) is (T|V^|5) = g[i B )2 x 
{(L\B z (v)\L) - (R\B z (v)\R) + (L\B z (v)\R)(R\L) 
-(R\B z (r)\L)(L\R)}, where the singlet state is \S) = 
[\L(1))\R(2)) + \L(2))\R(2))} (| U) - | |T»/2, the triplet 
state is \T) = [|L(l))|fl(2)> - \L(2))\R(2))] (| U) 
+ | 4T))/2, and the Zeeman coupling is Vz = 
g^B [B z {vx)Sx z + -B 2 (r 2 )52 Z ]. Here we have assumed 
that the external field is along the z direction. \L) and 
\R) are the left and right quantum dot ground states. 
Therefore, the inhomogeneous magnetic field couples the 
singlet and the S z = triplet states so that two of the 
eigenstates of the system are no longer the eigenstates 
of the total spin. When the overlap between the two 
ground state wavefunctions is small, the main contribu- 
tion to this coupling comes from the average field differ- 
ence between the two quantum dots. Below we explore 
the consequences of this loss of symmetry. 

When the overlap of the electronic wavefunctions is 
small, we can assume that the electron orbital degrees of 
freedom are frozen, so that an effective spin Hamiltonian 



quite faithfully describes the two-spin system 14 



H s = JSi • S 2 + -fiS lz + -f 2 S : 



2z : 



(2) 



where 71 = g^ B (L\B z (r)\L) and 72 = gfJ, B (R\B z (r)\R) 
are local Zeeman couplings due to applied or stray mag- 
netic fields. Here we have implicitly assumed, based on 
the small interdot wavefunction overlap, that the two 
spins are distinguishable, with spin 1 on the left dot and 
spin 2 on the right dot. We have also assumed that the 
field is entirely along the z direction. Whether such a 
choice is reasonable will be discussed later. 

Hamiltonian (|^) can be expressed in the basis of four 
two-spin states | TT), I T4), I IT), an d I 44), an( l its 
eigenstates can be easily obtained. The two polarized 
states are decoupled from the other two, which are mix- 
tures of singlet and S z — triplet states: \ipx) = | 4T), 
\>h) = I 44), 1^3) = Ci| T4) + Pa| 4T), and |^) = c 2 | T4 
) — ci I 4T)- Here the coefficients c\ and C2 satisfy the 
relations c 2 /ci = yj\ + (S/2J) 2 - S/2J and c\ + c\ = 1, 
where 5 — 71 — 72 represents the field inhomogeneity. 
The energies of the latter two states are also shifted from 
those of singlet and triplet states (with 6 being the shift): 
E x = J+A, E 2 = J- A, £ 3 = -J+V4J 2 + 6 2 = J+S, 
and Ex = -J-^AJ 2 + 5 2 = -3J-6, where A = 71+72 
represents the average magnetic field. An important 
question here is whether these mixtures and shifts will 
cause any error in quantum computation in the schemes 
based on the exchange interaction. After all, the swap ac- 
tion in these models depends on the perfect phase match- 
ing in the evolution of singlet and triplet states, as we 
will show below. Since swap operation is an essential 
component of a spin-based QDQC and several other ar- 
chitectures, we need to precisely quantify the effects of 
mixtures in singlet and triplet states on the swap action. 
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To determine whether swap is affected, we explore 
whether a product state of spins 1 and 2 will evolve into 
a product state (pure states for both spins) again. Our 
strategy here is to calculate the Schmidt number of the 
two-spin state. We will evolve our two-spin state, calcu- 
late its density matrix, trace out the second spin so that 
the density matrix of the first spin is left. We can then 
look for the eigenvalues of this density matrix. If at some- 
time it has only one non-vanishing eigenvalue, the state 
of the first spin is pure. We can then find out whether 
this pure state corresponds to a swapped state. 

Our initial state is a product state given by 



|0(O)) = (ai|T) + aaU))G9i|T)+AU)) 



(3) 



If the two electrons are located in two well-separated 
quantum dots in the beginning, the above product state 
does not violate the antisymmetry requirement of a 
two-fermion state. This state can be expanded in the 
basis of the eigenstates of Hamiltonian (|J). It then 
evolves under Hamiltonian (g). The two-spin state 
at time t takes the form \<f>(t)) = a x {5\e~ lEltjh \ TT 
> + a 2 Ae- iB »*/»| a) + ( ai fc Cl + a 2 {3ic 2 )e- iE * t / h \ U 
) + (ai/3 2 c 2 — a 2 PiCi)e~ lEit / h \ If). The corresponding 
density matrix for the first spin can then be calculated 
straightforwardly: Pl = Ti 2 {pi 2 } = Tr 2 {|<£(i)) (<£(*)!}• 
The matrix elements of p\ are all functions of the con- 
stants a%, a 2 , Pi, /3 2 , ci, c 2 , field inhomogeneity S, ex- 
change coupling J, and time t (pin is also a function of 
the average field A). The eigenvalue equation for p± is 

^ 2 - (pin + Pni) A + (pmPui - |pm I 2 ) = • ( 4 ) 

To have a pure state for the first spin, which means 
that only one eigenvalue of the density matrix p\ is non- 
vanishing, the last term on the left-hand-side of the above 
equation has to vanish: 



PmPm ~ \pimY 



o. 



(5) 



The explicit expression for Eq. (||) is quite compli- 
cated, so let us first look at the simple uniform field 
situation, when 8 = and c\ = c 2 = l/y/2. Equa- 
tion (||) then becomes \ai(3ia 2 p 2 — j[(ai(3 2 + a 2 (i\) + 
{a^-a^e^Ka^ + a^D-ia^-a^e^W 2 = 0, 
where 9 = AJt. The solutions here are e %e = ±1. When 
e w = 1, \Mt)} = ai| T> + a 2 e iAt \ |), the state of 
the first spin returns to its initial state with a phase 
shift between the two coefficients. When e lB = —1, 
\4>i (t)) = /3ie zAt | t) +(3 2 \ I), the swap is achieved with the 
exception of an additional phase that can be corrected 
easily with a single-spin operation. For example, a pulse 
sequence can be constructed based on the phase-shifted 
swap U psw to produce a CNOT gate 



Ucnot 



e"4 



cr 2y g»(f + - 



xUS 



■w 



Physically, a uniform field means that the Zeeman cou- 
pling couples to the total spin (including both electron 
spins), so that the Zeeman term commutes with the ex- 
change term in the Hamiltonian (||), and therefore does 
not change the eigenstates. The shifts in the energy lev- 
els of the polarized states cause additional phase shift, 
but can be corrected by applying an opposite magnetic 
field with the same pulse shape, magnitude, and length. 
In summary, an external uniform magnetic field does not 
qualitatively change the proposed QC algorithm. Logis- 
tically it makes the QC operation more difficult because 
of the necessary correction pulses. 

If the magnetic field is inhomogeneous, 5^0. Let 
us look at a simple situation when the initial state is 
l<X«)) = |TI>, i-e. ai = f3 2 = l and a 2 = (3 l = 0. Eq. §) 
now takes the simple form of 



\ Cl c 2 (l-e l ")(ci + cie lti )\ 2 = 



(7) 



where 9 = (4 J + 25) t. When c\ ^ c 2 , as is the case when 
5^0, the only solution of this equation is e 10 = 1, which 
corresponds to the "return" operation with an additional 
phase: \<j)\(t)) = ot\\ T) + a 2 e %At \ [). The condition for 
an exact swap does not exist anymore. To find the best 
approximation to a swap (so that the state of first spin 
is as close as possible to being "down" in the current 
special case) , we calculate the minima of the expression 
in Eq. (Q) and find that e 10 — —1 does still produce a 
minimum in |(1 — e l0 )(c^ + c 2 e l0 )\ 2 (which is, however, no 
longer zero). When this condition is satisfied, the density 
matrix of the first spin is 

= TTTaU)U I + 7^3 1 TXT I , (8) 



1 + x 1 



1 + 
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where x — 5/(2J). Therefore, the state of the first spin 
cannot be simultaneously pure and the same as the initial 
state of the second spin. The state of the first spin will 
remain mixed (when it is close to a swap) with the state 
of the second spin. According to Eq. (||) , there would be 
an error of the magnitude x 2 if we perform such a "swap" 
operation, which needs to be corrected. 

For GaAs-based QDQC, the Zeeman splitting is gps ~ 
2.55 x 10~ 2 meV/Tesla, while a single Bohr magneton 
produces a field of about 5 Gauss at 1 nm distance and 
almost nothing at a distance of 100 nm in GaAs. This 
difference in the magnetic field leads to a x 2 in the or- 
der of 10~ 6 , which is right within the capability of cur- 
rently available error correction schemes. Since the field 
involved in 8 is the average value over an entire quan- 
tum dot, local magnetic field inhomogeneity caused by 
impurities should not cause intractable errors. 

Mathematically, inhomogeneity in the magnetic field 
means that a part of the Hamiltonian (^|) does not com- 
mute with the exchange term JSi • S 2l thus the eigen- 
states will change, leading to errors in swap, which will 
cause errors in CNOT in the proposed sequential scheme 
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P). In the parallel pulse scheme |T(J, conditional phase 
shifts can be produced directly from the Hamiltonian (||) 
(with precisely controlled field inhomogeneity) , circum- 
venting the swap actions. However, stray fields, which 
may be invariably present in realistic solid state QDQC 
architectures due to magnetic impurities/imperfections, 
can still produce errors which need to be corrected (even 
in solely exchange-based models which have no applied 
magnetic field JIJ]). Furthermore, the impossibility of 
exact swap does make the transfer of spins through swap 
an error-prone process in the presence of an inhomoge- 
neous field. 

Finally, inhomogeneity in the external magnetic field 
is not only a possible error source in swap and other 
two-qubit operations, but also an important factor in 
the single-qubit operations determined solely by the Zee- 
man coupling. Here the main concern is the precise def- 
inition of the direction of the quantization axis for the 
electron spins. For example, if spin-up and -down along 
z direction correspond to the two states of our qubit, 
then the z axis is our quantization axis which provides 
a reference frame for all the single-qubit operations. In 
our calculations the magnetic field is purely along the 
z direction. The only spatial variation is in its mag- 
nitude. However, according to Maxwell equations, in 
a steady state, V x B = 0. Thus, if dB z /dx ^ 0, 
dB x /dz = dB z /dx 7^ 0. If we have a perfect two di- 
mensional QD system lying entirely in the x-y plane, we 
can assume that B x happens to vanish in the plane of 
the quantum dots. However, using the example of GaAs 
quantum well QD's, the thickness of the well is gener- 
ally around 5—10 nm, which characterizes the typical 
z-width of the QD's. Therefore, to have a vanishing Zee- 
man coupling along x direction, the B x field has to van- 
ish along the middle plane of the quantum well and be 
an odd function in the z (growth) direction, which is a 
quite stringent constraint. If we choose z-direction as our 
quantizing axis, a finite B x field, which may be unavoid- 
able, would tend to flip the spins. For a 1 Gauss B x 
field, this flipping rate is about 0.65 MHz, correspond- 
ing to a flipping time of 1.5 /is. Considering that the 
gates should be operated as slowly as possible to satisfy 
the adiabatic condition, this is a constraint one has to 
take into consideration. In fact, "decoherence" induced 
by a fluctuating inhomogeneous transverse field may turn 
out to be an important extrinsic decoherence channel for 
QDQC operations. 
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